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The study of the mean-field static solution of the Random Blume-Emery-Griffiths-Capel model, an 
Ising-spin lattice gas with quenched random magnetic interaction, is performed. The model exhibits 
a paramagnetic phase, described by a stable Replica Symmetric solution. When the temperature 
is decreased or the density increases, the system undergoes a phase transition to a Full Replica 
Symmetry Breaking spin-glass phase. The nature of the transition can be either of the second order 
(like in the Sherrington-Kirkpatrick model) or, at temperature below a given critical value, of the 
first order in the Ehrenfest sense, with a discontinuous jump of the order parameter and accompanied 
by a latent heat. In this last case coexistence of phases takes place. The thermodynamics is worked 
out in the Full Replica Symmetry Breaking scheme, and the relative Parisi equations are solved 
using a pseudo-spectral method down to zero temperature. 



Since its discovery, the spin glass (SG) phase has 
played and still plays a fundamental role in the inves- 
tigation and understanding of many basic properties of 
disordered and complex systems. The analysis of the 
mean-field approximation of theoretical models display- 
ing such a phase has revealed different possible scenarios, 
including different kinds of transition from the paramag- 
netic phase to the SG phase, as well as different kinds of 
SG phases. Most of the work, however, has been concen- 
trated on just two scenarios. 

In order of appearance in literature the first scenario is 
described by a Full Replica Symmetry Breaking (FRSB) 
solution characterized by a continuous order parame- 
ter function^ which continuously grows from zero by 
crossing the transition. The prototype model is the 
Sherrington-Kirkpatrick (SK) model, ^ a fully connected 
Ising-spin model with quenched random magnetic inter- 
actions. 

The second scenario, initially introduced by Derrida 
by means of the Random Energy Model (REM),'^ pro- 
vides a transition with a jump in the order parameter 
to a stable low temperature phase in which the replica 
symmetry is spontaneously broken only once. The order 
parameter is a step function taking two values gmin and 
the so-called Edwards- Anderson order parameter (Zea, 
(else said self -overlap), with q^i^ < gEA- In the paramag- 
netic phase they are both equal to zero. At the transition, 
gEA grows to a value larger than qmin- No discontinuity 
appear, however, in the thermodynamic functions. Actu- 
ally, at the transition to the one step Replica Symmetry 
Breaking (IRSB) SG phase, the Edwards- Anderson or- 
der parameter qea can either grow continuously from 
zero or jump discontinuously to a finite value. The 
first case of this second scenario includes Potts-glasses 
with three or four states;^ the spherical p-spin spin-glass 
model in strong magnetic field^ and some inhomogeneous 
spherical p-spin model with a mixture of p = 2 and 
p > 3 interactions The latter case includes, instead, 
Potts-glasses with more than four states^, quadrupolar 



glass models;^ p-spin interaction spin-glass models with 
p > 2iiSiii4iS and the spherical p-spin spin-glass model 
in weak magnetic field. ^ Because of the discontinuity of 
the overlap parameter across the transition the models 
belonging to this second case, often referred to as "dis- 
continuous spin glasses" (see, for instance, Ref. {IS]), 
have been widely investigated in the last years because 
of their relevance for the structural glass transition^iiii^ 
In all scenarios discussed above no latent heat occurs, 
i.e. the phase transition is continuous in the Ehrenfest 
sense. In this paper, instead, we consider a spin glass 
model undergoing (below a given critical point) a true 
thermodynamic first order phase transition between a 
paramagnetic (PM) and a Full Replica Symmetry Break- 
ing (FRSB) SG phase, presenting coexistence of phases 
and latent heat, completing the work presented in Ref. 

m 

Such a model is a generalization of the Blume-Emery- 
Griffiths-Capel (BEGC) modeliSiii for the A transition 
and the phase separation in the mixtures of He^-He^ in 
a crystal field, in which a quenched disordered magnetic 
interaction is introduced. In that case the phase diagram 
consisted of a fluid phase, a super- fluid one and a mixed 
phase (see Fig. Q] for a pictorial representation of the 
diagram). 

The first study of a spin-glass model undergoing a gen- 
uine first order thermodynamic transition is the Ghatak- 
Sherrington model (GS)fiS, a simplified version of the 
model under current investigation. In this spin- 1 model 
in crystal field, no biquadratic coupling is considered: the 
Replica Symmetric (RS) solution and its stability have 
been carried out by Lage and Almeida^^ Mottishaw and 
SherringtonSS and da Costa et almL, albeit not always 
with compatible results. There the evidence for a first 
order phase transition was found, in the neighborhood of 
the tricritical point. 

In the last few years some work has been devoted to 
the comprehension of the Random generalization of the 
BEGC model (RBEGC), mostly at the level of the RS 
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X = 1 - p 



FIG. 1: The phase diagram of the mean field model for He^- 
He'' mixture as originally studied by BEG. x is the density of 
He^ particles and p the density of He**. 



solution, that turns out to be unstable as soon as the 
transition takes place. Such an analysis has been per- 
formed in Refs. p^l23l | for what con cerns the lat tice gas 
version of the model, and in Refs. '24^ 25l26l27f for the 
spin-1 model. For sake of completeness we also report 
the study where, together with the random mag- 
netic interaction, also a random biquadratic interaction 
is considered. 

Quite recently, the spherical version of the present 
model has been worked out exhaustively for positive and 
null particle-particle interaction, both statically and dy- 
namically, by Caiazzo et. alw^t^. 

In section we present the model, also indicating 
the connection with different notation in the literature. 
In Sec. m the replica formalism for the model is re- 
called, the Replica Symmetric (RS) solution presented, 
together with its stability and low temperature analy- 
sis. In Sec. mil we adapt the variational method for the 
FRSB solution"*" to the model, we explain the numerical 
method of resolution applied and we discuss the practical 
necessity of using different gauges to express the thermo- 
dynamic functional in different parameters regions. In 
Sec. llVl we study the thermodynamic observables. From 
the behavior of the entropy versus temperature at fixed 
chemical potential we observe that below the tricritical 
point the transition involves latent heat. In Sec. IVI we 
display the phase diagrams in the parameters T, ^ and 
density. Sec. I VII contains our conclusions. 



I. THE RANDOM 
BLUME-EMERY-GRIFFITHS-CAPEL MODEL 

There exist two completely equivalent versions of this 
disordered model (at least as far as the statics is con- 
cerned). One version is a direct generalization of the 



BEGC model, with spin 1 variables (cTi = 1, 0, —1 on site 
i),^** the other one is a lattice gas of Ising spins (spin 
Si = 1,-1, with occupation numbers Ui = 0,1 on site 

In this paper we will use the second formulation 
described, in the mean-field approximation, by the 
HamiltonianS^ 



= - 2^ JijSiSjniTij-— y ^njUj-^y Ui-hy ^S^n^, 

i<j i<j i i 

(1) 

where the Ising spin glass lattice gas in an external mag- 
netic field is coupled to a spin reservoir by the chemical 
potential fi. K is the particle-particle coupling constant 
and h the external magnetic field. The magnetic interac- 
tion is described by quenched Gaussian random variables 
Jij, symmetric in i <-^- j, with mean Jy = Jq/N and vari- 
ance Jfj = Jjj^ — J^/N. Here and in the following the 
overline denotes average with respect to disorder. 

Just for completeness we also report the Hamiltonian 
in the original formulation: 



H = - ^ J^j(T^a■J - — ^(cri(Tj)2 + - ffi, 

i<j i<j i i 

(2) 

The transformation at = Si Ui, between the spin lattice 
gas dynamic variables and the spin-1 variables ai, and 
the transformation D — —fi — Tlog 2, between ^ and the 
crystal field D of the spin-1 system, allow for a perfect 
equivalence of two versions of the model. 

Some limiting cases of the model are the SK model^ 
(for ^/J — > oo) and the site frustrated percolation 
model^i (for K = —1 and J/1.1 00). 

In most cases studied up to now^^i^^i^Qi^^i^^i^'^i^'^i^^i^^i^'^ 
the analysis was mainly limited to the RS solution. The 
general picture which emerged was an instability of the 
RS solution below some transition line in the region of low 
temperature and large density of particles. The nature 
of the SG phase (IRSB of FRSB) and the type of tran- 
sition (SK-likei, p-spin-like^*^ or different!^) were, how- 
ever, not clear, also due to some "anomalous" properties 
at the RS level such as, for instance, complex stability 
eigenvalues. Despite this fact, the possibility of ap-spin- 
like transition has put new interest into this model as a 
possible, more realistic, model for the structural glasses, 
and its finite dimensions version has been numerically in- 
vestigated in a search for evidence of a structural glass 
transition scenario^SiS. Performing the quenched aver- 
ages in the most general Replica scheme of computation 
it has been possible to show that the stable solution in 
the mean-field case, is FRSB everywhere in the SG phase, 
where spins are frozen^ 
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II. THE REPLICA TRICK FOR THE 
THERMODYNAMICS OF DISORDERED 
SYSTEMS 

For any fixed (quenched) coupling realization J, the 
partition function of a system of N dynamic (annealed) 
variables a, is given h^^^^ 

Zn[J]^Tt, exp(-/?7^[J;a]) (3) 

and the quenched free energy per spin is 



J d[J] P[J] logZ^ 



[J] 

(4) 

where (■ • •) indicates the average over the couplings re- 
alizations. The thermodynamic limit of the free energy, 
— limiy^oo^og Zn[J] I N(i is well defined and equal to 
the quenched free energy / = limjv^oo /jv for almost any 
coupling realization J (self-average property). 

The analytic computation of the quenched free energy, 
i.e., of the average of the logarithm of the partition func- 
tion, is quite a difficult problem, even in simple cases as 
nearest neighbor one dimensional models. However, since 
the integer moments of the partition function are easier 
to compute, the standard method involves the so called 
"replica trick" by considering the annealed free energy 
j(n) of n non-interacting 'replicas' of the systenifSiS^iS^ 



N^oo N(3n 



log {ZN[j]y 



(5) 



The quenched free energy of the original system is then 
recovered as the continuation of f{n) down to the un- 
physical limit n = 0^ 



f 



lim lim 



{ZN[j]r-i 

N(3n 



lim/(n). (6) 

n— >0 



In the last equality we assumed that the replica limit and 
the thermodynamic limit can be exchanged. This proce- 
dure replaces the original interactions in the real space 
with couplings among different replicas. The interested 
reader can find a complete and detailed presentation of 
the replica method for disordered statistical mechanical 
systems in Refs. Q and [ssj . 

Applying the replica trick to the computation of the 
thermodynamic potential of the random BEGC model 
we find, in the saddle point approximation for large N, 



ZnW = exp {-nNpf [{p}, {m}, {q}]} 



4 n 



^Y^ql,-^l0gZ' 

n. n 



Z'^ exp{-/3i7'[{d},{m},{g}]} 
{S},w 



(7) 

(8) 
(9) 



n 

-PH' [{p}, {m}, {q}] = ^ (pkpa + (10) 



+ ^ SaUa if3Joma + (3h) 



2 1," 



qabSanaSbUb 



a=l 



where the sum in the one-site partition Z' is taken over all 
the possible values of the n spins Sa and the n occupation 
numbers Ha- Here and everywhere else in this paper we 
use the abbreviation 



K = K 



(11) 



The saddle point equations, coming from the extremiza- 
tion of Eq. ((SJ, give the self-consistent relations 



qab = {SanaSbUb) 

Pa = {no) 

nia = {SaUa) 



(12) 
(13) 
(14) 



where (...) is the average computed over the measure 
exp{-(3H'). 

The parameter pa represents the density of the 
replica a. In the thermodynamic limit this is equal to 
^/N^^rii — rii which, in the Hamiltonian is only 
coupled to the chemical potential fi, that is a replica in- 
dependent quantity. The density is therefore equal for 
each replica: pa ~ p, Va = 1, . . . , n. 

The same holds for rria which is coupled to the external 
field h: rua = m, Va = 1, . . . , n. It generally holds that 
one index quantities are replica invariant 



A. Replica Symmetric Solution and Stability of the 
Paramagnetic Phase 

The replica symmetric free energy is obtained by eval- 
uating Eq. ijHIl at qab = qo, for a ^ b, Pa = p and nia = m 
for every a, b and it reads^^ 



2 2 4 



q^~PJ dyPoiy) My) 



with 



eo 

Poiy) 
My) 



{p - qo) + (iKp + pn 



(3Kp + (3p - 



1 



1 



exp 



(/3J)"go 
2 

{v-m.h/J -h/jf 



2 go 



log(2 + 2 e'='° cosh (}Jy) 



(15) 
(16) 

(17) 
(18) 



The order parameters p, q^ and m satisfy the following 
self-consistency equations (see App. IVII for the details of 
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derivation) : 



P = 

10 



Po{y) p{y) 
dy ^0(2/) TO(y) 



with the foUowing definitions 

cosh /3Jy 



m(y) 



g-eo _|_ cosh/3Jy 

sinh /3 Jy 
g-eo _|_ cosh 13 Jy 



(19) 
(20) 
(21) 

(22) 
(23) 



The eigenvalues of the Hessian of Eq. computed in 
the RS approximation are derived in appendix I VII where 
the stabihty analysis is carried out. 

In the case Jo = h = only the paramagnetic solution 
ilo ~ 0) is stable, so that Eqs. ((T^ - (PT|l reduce to 



go -0 



1 



m = 



(24) 



and the eigenvalues of independent modes are 

Ao = iPJf [1 - ipJfp'] (25) 
Ai=P\K\[l-PK{l-p)p\ (26) 

where Aq is connected to Sq Sq fluctuations, whereas Ai 
is the stability eigenvalue of the density-density fluctua- 
tions. In the above formulas we have considered both the 
case in which > and K < 0. The last one occurring 
only if the particle-particle interaction K is negative and 
when T is bigger than —1/{2K). 

The lines Aq = and Ai = delimit the stability 
region for the RS solution on the phase diagram. In the 
T-p phase diagram the stable region is for 



T> Jp 

T>^p{l-p) 



corresponding, respectively, to Ao > and Ai > 0. 

For any value of K, there is one intersection between 
Aq = and Ai = 0, namely the tricritical point: 




J 

Pc 

J 



J 

2K 





3 


K 


(- 


2 


+ 7 + 


'K 








+ 


log ( — 
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\Pc 




where /ic is obtained from the paramagnetic expression 
Ea for p. In table d we list some values of interest of 



the tricritical values and in Fig. |21we plot their behavior 
as function of the particle-particle coupling constant. 



K/J 


00 


1 





-1 


— CX3 




1 


1/2 


1/3 


0.21922 





fJ.c/J 


—00 


-1 


-0.73105 


-0.55923 






TABLE I: Critical values of the thermodynamic parameters 
for some specific particle-particle interaction values. Param- 
eters are expressed in units of J. 




FIG. 2: The values of the thermodynamic parameters (in 
units of J) at the tricritical point as a function of the particle- 
particle coupling constant. The critical density is equal to Tc. 



As K decreases, the critical temperature goes to zero. 
As we will see in the following, in Sec. where we study 
the phase diagrams, this implies that the phase diagram 
region of phase coexistence progressively reduces itself. 
The critical value of the chemical potential grows to zero 
as K tends to —00: in order to contrast large particle 
repulsion a larger chemical potential (in the limit non- 
negative) is needed. 



1. Thermodynamic Observables in the Replica Symmetric 
Solution 



The internal energy and the entropy take the form: 



u = p — Jo m — n m — ^ p -\ — (-jl) 



m\ ,2 n 
s = -. — (p-qo) - p 6o 



(32) 

/oc 
dy Poiy) [My)-ym{y)] 
-00 

Integrating by part we can derive: 

dy Po{y) y m{y) ^ (3J{p - qo) (33) 
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The entropy, then, becomes 

-[p-qo) -pBo 



(34) 



4 

-{(iJfq^ {p - go) + PJ I Ay Poiy) Mv) 



B. Low Temperature Behavior of the RS Solution 

It can be of some help for later purposes and also to 
make a comparison with existing results, to probe how 
the observables behave in the very low temperature limit, 
even in this case for which the SG solution is unsta- 
ble. We simplify the discussion to the case Jq = h = 0. 
Putting z — y/y'qo as integration variable in the above 
expressions the entropy can be written as 



+!3J 



{p-qof - p ©0 



(35) 



dz 



_<Po{z) - V?o z m(z)] 



We define 



C = pjip - qo) 



-z'/2 



2tt 



(36) 
(37) 

{'/'o(z) - V9o z Th{z)} (38) 



We consider two cases separately in the zero tempera- 
ture limit. 



Full system: a > 



J 



lim 0=1 



lim C = 

lim a = —= + ■ 
A/ = Go + 0(T) 

»-2l + °<^) 



2. Partially diluted system: a < 



lim p = p < 1 

hm C = C 

__C K^+ji 
lim a = a — 1 

P^oo 2 J 



(39) 

(40) 

(41) 
(42) 
(43) 



(44) 
(45) 

(46) 



with p and C given by the zero temperature self- 
consistency equations: 



(47) 
(48) 



p = erfc(- — 



C = W — exp -== 

V TT p V V2 P 



Notice that as a ^ o ne ob tains from the above equa- 
tions: p = 1 and C = 
Finally one gets 

A/ = p Oo + (1 - p) log 2 - C a log 2 -I- 0(T)(49) 



C 



+ (1 - p) log2- a a log2 + 0(T) (50) 



For fixed K and J there is, thus, a limiting (negative) 
value of the chemical potential below which the lattice 
will not be full at zero temperature. For 



J 



1 



< - 



K 
~1 



(51) 



the density (and the overlap) is equal to p < 1. 

In any case the asymptotic value of C is finite, yielding 



(52) 



III. A VARIATIONAL METHOD FOR THE 
FULL REPLICA SYMMETRY BREAKING 
SOLUTION IN DISORDERED SYSTEMS 

To clarify which kind of transition takes place we im- 
prove the study of the static properties of the SG phase of 
the mean- field RBEGC, making use of the FRSB Parisi 
Ansatz. The first thing to notice is that the stable SG 
phase is always of FRSB type. The transition between 
the PM phase and the SG phase can be either of the SK- 
type or, below a given Tc, discontinuous with a jump in 
the entropy and hence a latent heat. Moreover for a cer- 
tain range of parameters, the two phases, PM and SG, 
coexist. For any parameter choice we find no evidence 
for a p-spin-like transition with discontinuous order pa- 
rameter. 

The aim is to evaluate the n ^ limit in Eq. © 
with the Ansatz that the structure of the matrix Q fol- 
lows a FRSB scheme. In order to be as general as possi- 
ble, we shall use the RSB scheme introduced for the SK 
model by de Dominicis, Gabay and Qrland^SSiSi which 
besides the Edwards- Anderson order parameter— also in- 
volves the anomaly to the linear response function, oth- 
erwise called Sompolinsky's anomaly.^® The more usual 
Parisi's RSB scheme is, if necessary, eventually recovered 
by a proper gauge fixing, once that the limit of infinite 
number of replica symmetry breakings has been taken 
[see below Eq. ((7^ ]. 

By applying the RSB scheme infinite times and in- 
troducing the two functions qix) (overlap function) and 
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A (a;) (anomaly function), < x < 1, the free energy 
functional, Eq. ||SJ), becomes: 



Pf = P — + ^Jq m 



(53) 



(3J 



dx qix) A{x) ~(3J Ay Poiy) 0(0, y) 



where Poiy) is defined as 



Poiy) ^ 



exp ■ 



iy-ih + Jo m)/J)' 
2q{0) 



(54) 



(56) 



and 0(0, y) is the solution, evaluated at a; = 0, of the 
Parisi parabolic equation^ 

k^, y) = -^0"(x, y) + ^0'(x, y)2 , (55) 
with the boundary condition at a; = 1 

0(1, y) = 01 (y) = (/3J)-i log (2 + 26®^ cosh/3Jy) , 

and 

ei^^[p-g(l)]+/% + ifp) (57) 

= /3ifp + /3M-g(l)^ . 

The overlap qix), the density of occupied sites p and 
the anomaly A(x) are the order parametersi^ We have 
used the standard notation and denoted derivatives with 
respect to x by a dot and derivatives with respect t o y 
by a prime. In this notation Sompolinsky's A' in Ref. l38l 
becomes our TA. 

The Parisi equation (|55|l can be included into a free 
energy variational functional via the Lagrange multiplier 
P{x,y) and the initial condition at a; = 1, Eq. H5t)|) . via 
P(l,y)i^ The free energy functional is then 



/3/v - PI 
+ pJ 



(58) 



dyPil^y) [0(1,2;)- 01 (y)] 



pJ I dx 

'0 



<iyPix,y) 0(a;,j/) 

<ii^):lf, X Mx).,, .2 

-^^0 i^iv) — ^0 i^iy) 



with Poiy) and 0i(?/) defined in Eqs. 

By such a construction /v is stationary with respect to 
variations oiPix,y), P(l,y), 0(a;,y), 0(0, y), g(x), Aq(a;) 
and deriving with respect to p. Variations of P(a:, y) and 
P(l, y) simply give back Eqs. (|S5|l and l(5H|) . Stationarity 
with respect to variations of 0(a:,?/) and 0(0,?/) leads to 
a partial differential equation for P(a:, y): 

^'(^,y) = ^^"(2:,y) + A(a;)[P(x,2/)0'(x,y)]' , (59) 



and to the boundary condition at a; = 
PiO.,y) ^ Poiy) . 



(60) 



Eventually, variations of /v with respect to qix), A(a;) 
and the derivative with respect to p lead to 



A(a;) = -pj [p - qil)] +/ dy P(a;, y) 0"(a;, y) (61) 

(62) 



p= / dyP(l,2/) 



dyPix,y) 0'(a;,y)2 
cosh/3 Jy 



e~®i + cosh /3Jy 



(63) 



with A(l) = 0, the anomaly at the shortest time-scale, 
corresponding to a; = 1, being zero by construction: the 
Fluctuation- Dissipation Theorem (FDT) holds at short 
time-scales. 

The Lagrange multiplier P(a;, y) represents the distri- 
bution of local fields. One may indeed associate a given 
overlap qix) with a time scale such that for times of 
order states with an overlap equal to qix) or greater 
can be reached by the system (these time-scales com- 
pletely decouple in the thermodynamic limit). In this 
picture the P(a;, y) becomes the probability distribution 
of frozen local fields y at the time scale labeled by xM> 

For the numerical treatment of the FRSB equations 
and also to allow for a clearer physical interpretation of 
the functions that we are analyzing, we define the local 
magnetization m(a;,y) = 0'(a:,y), whose differential an- 
tiparabolic equation we derive from Eqs. 1551) . 156(1 as 

rnix,y) = -|m"(a;,y) + A(a;) m(x, y) m'ix,y) , (64) 



m(l,y) = mi(y) 



sinh(/3y) 



e^Qi + cosh(/3y) 



(65) 



The average equilibrium magnetization to can, then, 
be computed in terms of the local magnetization at a; 0: 



a/ 

dh 



dyP(0,y) m(0,y) . (66) 



Another useful relation for the numerical evaluation of 
the order parameters and the thermodynamic functions 
built on them is got from the computation of the term 
da; J^^dy Pix,y)4>ix,y) in two different ways: once 
using Eq. (|^ and, in the other case, applying Eq. l(5(l 
and integrating by part. From the comparison of the two 
results the following equation is obtained: 



dy P(l,y) 0(1, y)- / dy P(0,y) 0(0, y) (67) 



1 

2 Jo 



dx qix) A(a;) 



Deriving Eq. H61|) with respect to x yields 

dy P(a;,y) m'ix,y)^ = 1 (68) 
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valid for every x and guaranteeing the marginal stability 
of the FRSB Ansatz. 

The coupled equations (|59|l . 1)64(1 . with border condi- 
tions H6U|) . Ht)5|l are the FRSB equations. In the following 
we are going to show how they can be numerically solved 
and what are the gauges we will use in order to get the 
most general SG phase and the order parameters charac- 
terizing the transition to it. 

Differentiating once more Eq. H68() one finds 

_ A _ /^dy Pix,y) m"{x,yf ^^^^ 
9 /^dj/ P(x,y) m"'(a;,?/)2 

that will become useful is the following when we will dis- 
cuss the choice of the gauge to perform the numerical 
computation (see Sec. IIII B|) . 



A. Numerical Integration of the FRSB Equations: 
the Pseudo-Spectral method 

In order to study the low temperature regime of the 
RBEGC in the limit of a large number of clauses we have 
numerically integrated the FRSB equations H59|l - Ht)4|) to 
determine q{x), P{x,y) and m{x,y). We followed the 
iterative scheme of Refs. ^40 41], but with an improved 
numerical method which allows for very accurate results 
for all temperatures (see Refs. [i^li^ 'l. 

We start from an initial guess for q{x), A(x) and p then 
m{x,y), P(x,y) and the associated q{x) are computed in 
the order as: 

1. Compute m[x,y) integrating from x = 1 to x = Q 
Eqs. H64|) with initial condition H()5I) . 

2. Compute P{x,y) integrating from a; = to a; = 1 
Eqs. lfCT|) with initial condition 

3. Compute q{x), A{x) and p using Eqs. Ht)2l()3|l . 

The steps 1 ^ 2 — > 3 are repeated until a reason- 
able convergence is reached, typically we require a mean 
square error on g, P and m of the order 0(10^®) and 
we checked that the identities H67|l. H81|l and 1 — T — q 
were satisfied to this precision as well. The number of 
iterations necessary are a few hundreds. The core of the 
integration scheme is the integration of the partial differ- 
ential equations (|59|) and Ht)4|) . In previous works this was 
carried out through direct integration in the real space 
which requires a large grid mesh to obtain precise re- 
sults. To overcome this problem we move to the Fourier 
space where we can apply a pseudo-spectral^ method of 
integration. 

Indicating by FT[o](a;, A:) the Fourier transform of 
function o{x, y), 

FT[o](x,/c) = — / dy e'^^y o{x,y) , (70) 



the FRSB Eq. (|56|l , written in terms of the wave number 
fc, becomes 

rn[x, k) = /c^iM m{x, k) + ik ^!^FT[m\x, k) (71) 
and the FRSB Eq. (jSHJ takes the form 

P{x,k) = -fc2^ P{x,k) + ik A{x) FT[P m\ix,k) 

(72) 

For each value of k these are ordinary differential equa- 
tions which can be integrated using standard methods. 
To avoid the time consuming calculation of the convolu- 
tion in the nonlinear terms we use the pseudo-spectral 
code on a grid mesh of x Ny points, covering the 
X interval [0,1] and the y interval [— j/max, ymax]- The 
truncation of the wave number may bring to anisotropic 
effects for large k. De-aliasing is, thus, performed by 
a Ny/2 truncation, which ensures a better isotropy of 
the numerical treatment. Eventually, the x integration 
is carried out by means of a third order Adam-Bashford 
scheme which reduces the number of fast Fourier trans- 
forms (FFT) callsi^ 

Typical values used are iV^, = 500 1000, Ny = 
1024-^-4096 and ymax = 24 -=-48. The difference between 
the values for and the values for Ny comes from the 
fact that, if the solution in y is smooth enough, only a few 
low wave numbers k are exited. The value of the param- 
eter j/max fixes the y range where the solution is assumed 
different from zero. Indeed, in the numerical algorithm is 
assumed that P{x,y) = m{x,y) = for |y| > y-max- This 
explains the rather large value of ^max used. 

B. Choice of Gauge 

The solution to the spin-glass mean-field models ob- 
tained using the scheme of Sommers is overconstrained 
and the functional expression of A(a;) can be, thus, cho- 
sen in different ways, selecting, in this way, a gauge for 
the order parameters. The usually studied gauges are 

A(a;) = —fiJ x q{x) Parisi gaugei (73) 

A(a;) = -A(0) Sommers gauged (74) 

where the anomaly is given by the stationarity Eq. I|61|) 
at x = 0. We rewrite it here in terms of the local mag- 
netization 

/oc 
dy P(0,y) m'(0,y) (75) 
-oo 

The anomaly at the largest time-scale is gauge invari- 
ant since it depends on the Edwards- Anderson parameter 
and on the derivative of the equilibrium local magnetiza- 
tion. A (a;) measures the violation of the linear response 
on the time-scale labeled by x. As such, it's a decreas- 
ing function of x, being zero at the shortest time-scale 
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(x = 1) where the dynamical relaxation of the system 
has not yet fallen out of equilibrium, and maximum at 
the longest time-scale x — 0. 

Other "intermediate" gauges can be introduced for nu- 
merical purposes, e.g. 



to be derived as 



A(a;) = -7A(0)(l-x)'^- 



7= 1,2,. 



(76) 



that, for 7 = 1, is the Sommers gauge. 

In practice, the most common choice is the Parisi gauge 
A(a;) — —(3 J x q{x), useful from a numerical point of 
view, since for x larger than a certain critical value Xc, 
at which the overlap function displays a cusp and reaches 
the plateau value qEA, tends to zero and the integra- 
tion domain can, thus, be reduced. For T — > 0, however, 
the overlap function in this gauge tends to a step func- 
tion, yielding, for some x, a diverging factor q{x) in both 
terms of the rhs of equations (|71|l and H72|l . 

For this reason the alternative Sommers gauge can be 
adopted A(a;) = — A(0), so that the non linear terms of 
the rhs of Eqs. (|71|l and (|72|) is automatically kept finite. 
The function q{x) comes out to be a smoother function 
of X in this gauge, not only continuous and monotonous, 
but also without any cusp at any x. When the external 
magnetic field is zero q is still divergent, but only exactly 
at cc = 0, thus yielding no sharp change in convexity, 
whereas for h ^ even this divergence disappears and the 
FRSB equations do not display any integration problem 
because of q{x) behavior. 



dq 



(78) 



2. Parisi gauge as fluctuation- dissipation ratio 

Following a dynamical interpretation of the order pa- 
rameters q{x) and A(x), respectively as spin-spin correla- 
tion function and anomaly in the susceptibility at time- 
scale X, the Parisi gauge is actually a rewriting of the 
Fluctuation-Dissipation Relation. That is the generaliza- 
tion of the FDT for aging systems dynamically stuck out 
of equilibrium, such as spin-glasses, where the coefficient 
between correlation and response functions is not T but 
some function of the time-scale on which the system is 
relaxing (often referred to as effective temperature). We 
saw in Eq. H89|l that A(a;) = xi^) ~ x(l) is the anomaly 
of the susceptibility with respect to the linear response 
value (i.e. the ZFC susceptibility). Thus, A(a;) is the 
anomalous response function, whereas q is the derivative 
of the correlation function and Tgg = 1/ {f3x) the effective 
temperature. The time dependence of Tcff is expressed 
through the time-scale index x. 



IV. THERMODYNAMIC OBSERVABLES 



1. Gauge dependent parameters and physical observables 

The Parisi gauge is, thus, advantageous down to tem- 
peratures where it starts bringing numerical instabilities. 
For lower temperatures the Sommers gauge (or some sim- 
ilar, e.g. Eq. H76|l ) stabilizes the results. To choose a 
gauge means to select the x behavior of the order param- 
eter functions will be different. This will not change, any- 
way, the physical quantities, such as density, Edwards- 
Anderson parameter (qea = 9(1)), energy, entropy, etc., 
that we are going to show in the next section, 

Also the overlap probability distribution P{q) is gauge 
invariant. However, one has to consider the fact that the 
definition P{q) = dx{q)/dq only holds in the Parisi gauge. 
To build a general definition of P{q) as the derivative 
of a cumulative distribution we can use Eq. H69() and 
introduce the cumulative function ^ 



TA{x) 
q{x) 



(77) 



where x here is the gauge-dependent RSB parameter and 
A and q are the gauge-dependent order parameter. In 
the Parisi gauge is ^(a;) = x, and we get back the usual 
definition. In the Sommers gauge is £,{x) — TA{0)/q{x) 
and ^(q) can be obtained, e.g. parametrically, using Eq. 
(|62|l . Thus, the expression of the overlap distribution has 



All thermodynamic quantities can be written in terms 
of the order parameter functions derived by solving Eqs. 
(|56|l . H59|) . namely overlap, density, anomaly, local mag- 
netization and distribution of local fields. 



3. Internal energy 

The internal energy u = d[3f /d[3 can be computed 
either taking the derivative of Eq. (and then breaking 
the replica symmetry infinite times) or directly deriving 
Eq. 

In the first case the energy comes out to be 



u = lim lim 



K + (3,P 1 ^ 2 r 1 2 



1 \ ^ 1 \ ^ /3 1 \ ^ 2 

-/l- > TOa-/i- VPa 7, la. 

T) ^ — ^ n ^ — ^ 2 n — 



K + PJ^ 2 Jo 2 , 
p TO — ti m — II p 



(79) 



+ ^<?(1)' + ^ f dx q{x) A(x) 



where Nb is the number of replica symmetry breakings. 
Otherwise, deriving Eq. (|53|) and using Eqs. H56lt)7|) . one 
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finds 



u = p - ^lp + —m 



(80) 



-J 



d2/f'(l,y) V m{l,y) 



The comparison between Eq. H79|) and Eq. H8U|) yields 
the relation 

2 [ dx q{x) A{x)=pJ q{l) [p~q{l)] (81) 
Jo 



Jo 2 , ^ 

m H m — 

J J 



dyP{l,y) y m{l,y) 



Taking the limit p — > 1, with K = fj, = and choosing 
the gauge A — —fiJ x q{x) (corresponding to the Parisi 
scheme for RSB) we get the SK formulas^. 



4- Entropy density 




FIG. 3: Compressibility versus temperature for K/J = 1, in 
units of J. The behavior at ^ above /ic/J = — 1 is smooth at 
the transition (right dotted line). At p.^ a discontinuity starts 
to develop (middle dashed curve) leading to a cusp for p < pc 
(left solid curve). 



The entropy density s — (3'^df/df3 can be expressed 
either as 

S = -pe,~^-^[p^qilf (82) 

/oo 
dy P{l,y) [cl){l,y)-ym{l,y)] 
-OO 

or, exploiting Eq. JSTJ, as 

^ - -^[p-g(l)f-pei-(/3J)2g(l) [p-q{l)] 



-(3 Jo - (3hm+ 2(3 J / dx q{x) A{x) 

Jo 

/•oo 

+/3J / dyP(l,y) 0(l,y) 



(83) 




0.4 0.5 0.6 0.7 



5. Compressibility 

The compressibility in terms of the density of occupied 
sites and its conjugated field, the chemical potential, can 
be expressed in the FRSB formulation as: 



_ 1 dp _ (3 
p^ dp p^ 

where we define 

pi(y) 



dy P{l,y) pliy) 



cosh{/3Jy) 



cosh{P Jy) 



(84) 



(85) 



In Fig. II V 51 we show the compressibility behavior versus 
temperature of the system with particle-particle interac- 
tion constant K = J for three values of the chemical 



FIG. 4: Field cooled and zero-field cooled susceptibility for 
K = 1. Above Pc — —1 FC (xo) and ZFC (xi) susceptibilities 
continuously reach the same value at the second order phase 
transition (xo = Xi = 1 in absence of external magnetic field) , 
whereas for p < pc both a discontinuity takes place and the 
value of xi never reaches the one of xo at the, now first order, 
phase transition. 



potential. Respectively above, at and below the critical 
value Pc below which the transition happens to be first 
order in the Ehrenfest sense. For p < p^ ^ cusp shows 
up and, decreasing the temperature, n shrinks to a much 
lower value as the transition point is crossed. 
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6. Susceptibility 

The magnetic susceptibility of the RBEGC model at 
equilibrium can be written either as 



Xo = PJ 



dx q{x) 



or, with the help of Eq. (|66|l . as 

Xo = ^= / dyP(0,y)m'(0,y) , 



(86) 



(87) 



whereas the susceptibility obtained when the system is 
constrained to stay in a single minimum of the free energy 
function comes out to be 



Xi=f3J [p~q{l)] 



(88) 



The equilibrium susceptibility is a function of the aver- 
age q = dx q{x) of the overlap over all possible values 
it can take at all time scales (0 < a; < 1). It corresponds 
to the Field-Cooled susceptibility. The second suscepti- 
bility, instead, only depends on the Edwards- Anderson 
parameter ^ea — 'z(l)- K physically expresses the self- 
overlap of configurations belonging to the same state. 
Indeed, at time scale n, i.e., the shortest time scale, the 
system has not yet visited but one metastable state, thus 
the response to a field perturbation only depends on the 
self-overlap. The experimental analogue of xi is the Zero 
Field Cooled susceptibility. As a matter of fact, also in 
that case the system remains in one single state during 
the cooling down to the SG phase, since it is not driven 
by any external field. 

From a dynamical point of view, the equilibrium sus- 
ceptibility Xo can otherwise be expressed in terms of 
the function A(a;) defined by Sompolinsky to encode the 
anomalous response to a field perturbation at large time 
scales (a; < 1, Xj; > Ti in the parametric representa- 
tion used so far). The anomaly function is a direct way 
to measure ergodicity breaking occurring in spin glasses 
(see Ref. i34.35,') . even at infinite time (x = 0). Defining 
the susceptibility function at the time-scale labeled by x 
as 



X{x) = pj[p - q{l)] + A(.x) 



(89) 



where A{x) is given by the stationarity Eq. (|(jl|l . the 
equilibrium susceptibility, Eq. H86I) . can be rewritten as 



XO - X(0) = fiJ[P - 9(1)] + A(0) 



(90) 



The anomaly A(0) can, thus, be interpreted as the dif- 
ference between the theoretical descriptions of the zero- 
field-cooled and the field-cooled susceptibility. 



A(0) - X(l) - X(0) - /5 



9(1)- 
PiqEA 



dx q{x) 



(91) 
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FIG. 5: Entropy density as a function of temperature for 
K = J. For fi < = — J the entropy is discontinuous at the 
transition temperature: a latent heat is produced/employed 
at the transition. 
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FIG. 6: p — q{l) as a function of temperature for K — J. 
For fi < jj,c = —J there is a discontinuity at the transition 
temperature. 



7. Free Energy Integral dy P{0, y) 0(0, y) 

The often occurring integral dy P(0, y) (f>{0, y) 
[see Eqs. and (|^ ] can be expressed as a function 

of the local magnetization m(0,y), numerically found as 
solution of coupled Eqs. H59|) . lt)4|) a; = 0. From the 
magnetization definition we write the identity 



(/.(0,y) = 0(0,0)+ / dy' m{0,y') , (92) 

that is valid for every y. Combining this with the identity 
(j){^,y) — 1^(0,?/) = Jq dx (t){x^y) and taking the limit 
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FIG. 7: p — p as a function of temperature for K — J. 
For fi < = —J J) < 1 and there is a discontinuity at the 
transition temperature. 



\y\ 3> 1, we get for the constant 0(0, 0) the value 

0(0,0)= hm |t Oi + lyl-i/ dxA{x) (93) 
\y\^oo I z Jo 

y 

dy' m(0,y') 



= rei-i fdxA{x)+ [ dy [l-m(0,y)] . 
^ Jo Jo 

This leads to the more convenient expression, for a nu- 
merical computation, 



dy P{0,y) 0(0, y) 



(94) 



= Tei-- dxA{x)+ / dy [l-m(0,y)] 
^ Jo Jo 



dy P(0,y) / dy' m(0,y') 

30 Jo 



1. Full lattice at zero temperature: a > 

For T ^ the last term of the entropy in Eq. 
goes as 

dy P(l,y)[0(l,y)-y m(l,y)] (95) 

-OO 

= ei(o) + o(e-i/r) 

where the boundary functions 0(1, y) and rn(l,y) are 
given by Eqs. The Edwards- Anderson parame- 

ter and the density both tend to one. 

At zero temperature there is one single ground state 
and the entropy is, therefore, zero. So that, from Eq. 
we get 



i^[p-q{l)f + {l^p)Q,\ (96) 



= lim 



Using Eq. 1)95(1 and equating the entropy expression 
at zero temperature to zero we see that q{l) and p go to 
their T — value in such a way that p ~ q{l) ^ T'' and 
l-p~r"',with s,s' > 1 

To analytically determine the value of s and s' is turns 
out to be quite complicated. However, from our numeri- 
cal computation, we find that, at low temperatures both 
the entropy and p~ q{\) reach the zero value as T^, thus 
s = 2. In order to satisfy formulas H82|) and H99|) . then, 
1~ p has to decrease as [s' = 3) , also consistent with 
our numerical data. 



2. Partially empty lattice: a <Q. 

The parameter Qi tends to —oo at zero temperature 
and q{l) and p are less than one, namely 

/•oo 

lim q{l) = lim p = p = 2 / dy P(l, y) (97) 



/3-»oo 

= 2 



,3^00 



-pK/J-^t/J 



dy Pil,y-pK/J-p/J) 







If 81 < 0, in general, e®^ cosh /3Jy will not be much 
larger than one for any value of y in the argument of 
0(1, y) and in y m(l,y). The zero T limit depends on 
the sign of a\y\ = \y\ + a/ J: 



A. Low Temperature Behavior in the FRSB SG 
Phase 



From the self-consistency Eqs. (|62I63|I we can get the 
T = limits of g(l) and p. As in the RS case previously 
reported (see Sec. Ill B|) . also for the stable solution that 
limit crucially depends from the sign of 0i at zero tem- 
perature. Once again we define a = Qi{T)/{j3J) and we 
sketch the results in the two qualitatively different cases: 



lim 0(1, y) = 



if a\y\ > 



(98) 



^log2 ifa|y|<0 



The analogue of Eq. H95|l is, in this case. 



13 J I dy P(l,y)[0(l,y)-ym(l,y)] (99) 
= (1-p) log2 + pei(0) + O(e-i/^) 
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Since the zero temperature density is less than one, there 
will be a fraction 1 — p of spins whose orientation is ir- 
relevant for measuring any observable. This brings to 
a degeneracy in the ground state, of 2^'^^"^^ equivalent 
configurations. Then, for the entropy density, we get 

s = lim \-i^[p-qil)f (100) 

+(p-p)ei + (i-p)iog2} 

= (l-p)log2. 

Also in this case we can see from numerical data at low 
temperature that the behavior of p — g(l) and s is T^, 
whereas p — p ~ T'^ is still consistent with our numerical 
datas. 

At zero temperature the free and internal energy are: 

u = f = -^p- P+ f dxq(x)A{x) (101) 
^ Jo 

where p is either 1 or p depending on the sign of Qi. 
V. PHASE DIAGRAMS 
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FIG. 8: T - p phase diagram of the RBEGC for K/J = 1. 
The dot marks the tricritical point pc/J ~ —1, Tc = J/2, 
Pc = 1/2. In the upper-left region we have the paramagnetic 
(PM) phase. In the upper-right region the FRSB spin glass 
(SG) phase. Following the iso-potential line at p/ J = —1.05 
one sees the jump in density at the first order phase transition. 



Analyzing the stability of the RS solution [see Eqs. 
()25I26|| ] one gets the critical lines 

1 - ipJp)^ = , (102) 
I- (3K{1- p)p = , (103) 

above which the only solution is the PM solution q{x) = 
for X e [0, 1], p = 1/[1 + e-®i] [see Eq. (EH)], stable for 
any value of AT. In the T—p plane, these are, respectively, 
the straight line and the left branch of the spinodal line 
shown in Fig. |S1 for A" = J. The two lines meet at the 
tricritical point lf^ - (pn|l . Between them, under the two 
hard line curves branching out from the tricritical point, 
there is a region of coexistence of phases (as indicated 
in the plot of Fig. The broken line curves are the 
first order transition lines. When p < pc, iu cooling at 
fixed chemical potential, the density of the system jumps 
discontinuously from the PM to the higher SG value. As 
an example the line at p = —1.05 is plotted. 

By crossing the critical line Eq. (|102|l above the tricrit- 
ical point (p > Pc, r > Tc, p > Pc) the system undergoes 
a continuous phase transition of the SK-type to a FRSB 
SG phase, with a non-trivial continuous order parameter 
function q{x) which smoothly grows from zero. 

Below the tricritical point the scenario is completely 
different, displaying a discontinuous transition from the 
PM phase to a FRSB SG phase with q{x) which dis- 
continuously jumps from zero to a non-trivial (continu- 
ous) function. At the critical temperature the entropy 
is discontinuous, see Fig. |S1 and hence a latent heat is 
involved in the transformation from the PM to the SG 
phase, implying that the transition is of the first order 
in the thermodynamic sense. The transition line is deter- 
mined by the free energy balance between the PM and 
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FIG. 9: p-p phase diagram of the RBEGC for K = 1. Above 
the tricritical point (star) the transition is second order. Two 
isothermal lines are shown for temperature T — 0.5, 0.6, i.e. 
at and above. 

the SG phase, and is shown as a broken line in the phase 
diagrams. The line yielded by Eq. (|103|) where the PM 
solution becomes unstable, and the equivalent line from 
the SG side are the spinodal lines. 

Also interesting is the p — p phase diagram represented 
in Figs. El El We indeed sec that the isothermal lines 
cross the instability lines with vanishing derivative (Fig. 
lion and hence a diverging compressibility k occurs cross- 
ing these lines. 

It can be shown that the first order transition line 
can be determined in the p — p phase diagram from the 
isothermal and spinodal lines by using a Maxwell con- 
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FIG. 10: /i-p phase diagram of the RBEGC for K = 1. From 
the tricritical point (star) two solid hnes come out: the spin- 
odal lines at which the paramagnetic (upper-left curve) and 
the spin glass (lower-right curve) phases cease to exist, even 
as metastable. The first order transition lines, also branching 
out of the tricritical point, are plotted as dashed curves. An 
isothermal line is shown, for a temperature T — 0.3 below 
Tc = 1/2. Along such curve a first order phase transition 
occurs. In the plot, also the metastable branches are shown, 
both in the RS PM phase and in the FRSB SG phase (broken 
curves continuing the full curves). They reach the spinodal 
lines with zero derivative. In this plane of conjugated ther- 
modynamic variables a Maxwell construction can be explicitly 
performed. 
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FIG. 11: T - ^ phase diagram of the RBEGC for K ^ 1. 
For low jj, (lower than -1.2 in this specific case) the code at 
T = - in the Sommers gauge - is very unstable. Therefore 
the zero-T transition points are obtained by fit. 

struction. In the region between the first order transition 
line and the spinodal hne the pure phase is metastable. 



Below the spinodal lines (in the T — p plane) no pure 
phase can exist and the system is in a mixture of PM 
and SG phase [phase coexistence). 

Eventually, the phase diagram in the T — /i plane, for 
K — \, IS shown in Fig. 1111 Since our code, even in 
'most suitable' gauge (see Sec. IIII Bll . is unstable in re- 
gion of both low temperature and low chemical potential 
we were not able to reliably compute the points belong- 
ing to the SG spinodal line and the first order phase 
transition line for T < 0.1 = Tc/5. The prolongation 
of the lines down to zero temperature are, hence, com- 
puted by fit. Therefore the zero-T transition are esti- 
mates rougher than the others. The first order line goes 
to /ij^' = —1.256 ± 0.009 and the spinodal one reaches 
zero at /Iq*^ — —1.465 ± 0.008. Nevertheless, it seems 
that we can rule out a reentrance to the PM phase as 

By varying K the scenario remains qualitatively un- 
changed [look at Ref J^. for the phase diagrams of the 
Ghatak-Sherrington modeliS [K = 0) and the frustrated 
Ising lattice gas model^2. [K — — J)]- The only effect of a 
strong repulsive particle-particle interaction is to increase 
the phase diagram zone where the empty system (p = 0) 
is the stable solution. In order to find further phases, 
e.g. an antiquadrupolar phasef23> a generalization of the 
present analysis to a two component magnetic modeljS 
including quenched disorder, has to be carried outsS 

If we now look at the qualitative reproduction of the 
phase diagram of the original BEG model, with ferromag- 
netic interaction, in Fig. ^one can imagine a correspon- 
dence between the PM phase of the disordered magnetic 
material and the fluid phase of the He'^-He'* mixture, and 
between the SG phase and the superfluid phase. The den- 
sity p corresponds to the He'' density and the density x 
of He^ particles is 1 — p. The second order transition 
line corresponds to the A transition points from fluid to 
superfluid. 



A. Transition Lines 

The starting point to analytically determine the tran- 
sition lines around the tricritical point is the expan- 
sion of Eq. JHJ) for small q and densities next to pc 
that we report in appendix IVII [Eq. p65(l ]. There the 
various coefficients are expressed in terms of the func- 
tion po = (1 + exp{—f}p — f}KT/J))~^ coinciding with 
the paramagnetic density of the system evaluated along 
the second order transition line Aq = (along which is 
P = Po = T/ J). The expansion parameters are the ele- 
ments of the overlap matrix qat and r = p — po (replica 
independent) . 

One can obtain the approximated analytical expression 
of the spinodal line at which the FRSB solution disap- 
pears and the first order transition line, at which the 
paramagnetic free energy overcomes the approximated 
spin glass free energy, around the tricritical point. 

In figure [T^ we show both the behavior of the spinodal 
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qi = -hT+]^^Je-5T^ (106) 

r = -^ + iv/7^«^ (107) 

Since 5T is always negative the above parameters are 
always positive, as far as they exist (e > ST"^). The 
solution breaks down at the spinodal line, therefore given 
by e = 5T2, i.e. 

11 = ^i,-T, + T +{T - Te)2 (108) 

From the comparison of the free energy values for the PM 
and the SG phases the approximated first order transition 
line turns out to be 

M = Mc-T, + T+ 11. 19(T-Te)2 (109) 



FIG. 12: The T-^ phase diagram around the tricritical point. 
The points displayed are those computed by directly solving 
the FRSB equations at given T and probing in where the 
SG solution disappears (SG spinodal) and when the SG free 
energy becomes lower than the PM one (first order transition). 
The lines are obtained, instead, from the free energy expanded 
to the second order in 5T. Also the second order transition 
and the PM spinodal line are displayed. 



line of the spin glass phase and the behavior of the first 
order transition line for the case K = J , va the neigh- 
borhood of the tricritical point in the T-^ plane. The 
points in the plot are obtained by numerically solving 
the FRSB differential equations. Notice that, at the or- 
der of approximation used, the first order line displays a 
reentrance that does not occur, instead, performing the 
exact computation. 

Notice that, at the present degree of approximation, 
already at such small distance from (/ic, Tc) the first order 
line displays a reentrance, whereas the exact computation 
shows that this does not take place. 

In the following we consider /i/J — > /i and T/ J — > T 
and we express everything in terms of the small quantities 
5iJL = ji — He ST = H — Tc and 



_ 9mo(T) 



(104) 



This last auxiliary variable represents the distance from 
the second order phase transition line. If e < 0, in the T- 
/i plane, we are above the second order line (PM phase), 
otherwise below (SG phase). The function iio{T) is de- 
fined as: 



1 



(105) 



The expressions of r 
of q{x), are 



P- Pq 



and qi, the highest value 



VI. CONCLUSIONS 

In the present paper we have shown, in some detail, the 
properties of the Random Blumc-Emcry-Griffiths-Capel 
model. Such a model can be seen both as a 1-spin model 
or a Ising-spin model on a lattice gas. The two formula- 
tions are completely equivalent, at least from a static 
point of view. In the second representation, that we 
adopted, the interactions involved are a quenched ran- 
dom magnetic coupling Jy between spins and an at- 
tractive/repulsive coupling between (full) sites (K). The 
system is embedded in a reservoir and the exchange of 
particles is controlled by the chemical potential ii. The 
quenched random interaction is the source of a spin-glass 
phase at low temperature, provided that the chemical 
potential is large enough. If, on the contrary, /i is lower 
than a certain value, the system happen to be always in 
the paramagnetic phase, becoming progressively empty 
as T — > 0. We analytically studied the system in the 
mean field approximation. As already shown by the au- 
thors in Ref. the qualitative features of the system 
do not depend of the value of K, nor on the coupling 
being repulsive, attractive or zero. 

The external parameters are the temperature and the 
chemical potential. In a certain region of the phase dia- 
gram (see figure in Sec. for the T — fi diagram at 
particle-particle interaction K — J, or, else, figures |H1 El 
and llO|l the system undergoes a second order transition 
varying T or /i, down to a given 'tricritical' point at which 
the second order line ends. Indeed, the main feature of 
this model is that for low temperature, or high chem- 
ical potential, a continuous transition (in the Ehrenfest 
sense) from a pure PM phase to a pure SG phase does not 
occur anymore. In its place a first order transition takes 
place, with consumption/production of latent heat (see 
Fig. Oin Sec. lIVp and the appearance of a region of the 
parameter space where the two phases (PM and SG) do 
coexist. The SG phase comes out to be stable exclusively 
in the FRSB scheme of computation. We can, thus, rule 
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out the existence of a glass-like phase, by this meaning a 
phase with one step RSB in the static, corresponding to 
a dynamic decoupling of the characteristic time scales of 
the processes involved in two time sectors. Only one kind 
of spin glass phase exist in the frozen phase and this is the 
one typical of mean-field models for amorphous systems 
(i.e. spin-glasses in the proper sense). In order to recover 
a system with a stable IRSB phase and displaying a first 
order transition a lattice gas of spins interacting through 
multiple spins interaction (i.e. a p-spin spin glass on lat- 
tice gas) has to be considered (see e.g. Ref. 48]). It is 
interesting to notice that a discontinuous transition be- 
tween liquid and glass, with coexistence of phases, has 
been recently found in lattice heteropolymers with ran- 
dom interactions^Si^. We also mention that generalizing 
the present model to a two component magnetic modelS, 
and including quenched disorder, further phases can be 
found, e.g. an antiquadrupolar phase. 

Computing the state of the system down to very low 
temperature (including zero for not extremely low /i val- 
ues) it has been possible to see that the reentrance dis- 
played in the T-fi phase diagram in the Replica Symmet- 
ric approximation is just an artifact (see e.g. |23|, but 
also the 'small g' expansion in Appendix B). There is 
no /^-range for which lowering the temperature the spin 
glass can transform itself in a paramagnet. 

We have discussed in detail the numerical method that 
we use to solve the antiparabolic differential Parisi equa- 
tion for the present model, allowing us to compute the 
overlap order parameter and all the thermodynamic ob- 
servables in the SG phase. In particular we face the prob- 
lem of making the code converge at very low temperature, 
including zero. Section imi is dedicated to the presenta- 
tion and explanation of the variational approach to the 
problem of computing the FRSB anti-parabolic differen- 
tial equations and the pseudo-spectral method employed 
to solve them. 

Besides the numerical resolution of the FRSB equa- 
tions we have also studied the phase diagrams around 
the tricritical point making an expansion for small over- 
lap values and densities next to the density at the tricriti- 
cal point (see Sec. IV Al and aDDendix lVI|) . The expansion 
formalism is valid all along, and around, the second order 
transition line. The only interesting point around which 
is worth probing the parameter space is, however, the 
tricritical point. From there one can build analytical ap- 
proximated expressions for the spinodal line of the spin 
glass phase and for the first order transition line as well. 
Above the tricritical point the transition is always second 
order and qualitatively identical to the phase transition 
taking place in the Sherrington-Kirkpatrick model, that 
is easily recovered as limit of the present model (density 
one, zero chemical potential, full lattice). 



APPENDIX A: Stability Analysis in the Replica 
Formalism 

Rescaling /3J ^ /?, Jq/J J, K/J K, h/J -> h, 
fi/J —^ /X and, consequently K K/J we write the 
replica thermodynamic potential ||SJ) as 

G({p}, {m}, {q}) ^ n/3/({p}, {m}, {q}) (110) 

\ " 2 , \ " 2 , /^^ \ " 2 1 ryl 

a—1 a—1 a^b 

with Z' and H' given respectively in Eqs. lO, H1U|) and 
repeated here for clarity: 

Z'^ Yl exp{-(3H'[{d},{m},{q}]} (111) 

{S},{n} 

-fiH' [{p}, {m}, {q}] = Y,n, [pKpa + (112) 

a 

-I- ^ SaUa iPJoma + 0h) + — ^ qabSanaSbUb 
o=l a^b 

The variation of G with respect to the replica param- 
eters is 

n 

6G{{p},{m},{q}) = pKY,PaSpa (113) 

a=l 

ma 5ma + -^'^Qab Sqab - ^ I3K {Ua) 5pa 
a—1 a^b a—1 

+ ^ PJo {Sana) Snia + "^51 (SanaSbUb) Sqab 
a—1 a^b 

where the average (. . .) is performed with the measure 
given by Eq. (|112|) . From SG the saddle point equations 
are immediately derived: 

qab = (SanaSbUb) (114) 
Pa = (Ua) (115) 
rria = (SaTIa) (116) 

Eventually the fluctuations functional is equal to 

n 

S'G{{p},{m},{q})^PKY.(^P-^' ("7) 

a=l 

" ^32 

a—1 a^b 
n n 

-Y,(3k5{{na)) Spa + Y,I^JoS{{Sana}) Sma 

a—1 a—1 

+ -;^'YS{{SanaSbnb)) Sqab, 

a^b 
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where 

5{0a^...0a^) =5 — (118) 

and 

n 

5{-l3H')= l3kY,n,5pa (119) 

a=l 

a 

a=l a^b 

In H117|l we have, thus, six kinds of terms: the 
quadratic ones and the mixed ones whose coefRcients are 
given in table ITU 

The functional S^G, that has to be definite positive 
in order for the solution around which the expansion is 
performed to be stable, takes the form 

d'G = P'Y.^Sqab)'+PKY.^Spa)' (120) 

ab a 

+PJo^{Smaf + ^ Sqab A(^ab){cd) Sqcd 
a {ab),(cd) 

+ ^ Sqab D(ab)c Spc + ^ Sqab E(ab)c Suic 
{ab),c (afc).c 
'. I 





Second order expansion term S'^G 


fluctuation term 


coefiicient 


Sqab Sqcd 


A(ab)(cd) ~ —P^ {{naSanbSbncScTldSd) — {UaSaTlbSb) ijlcS cTtdS d)) 


Spa Spc 


Bac = -iPKf {{naUc) - (Ua) {Uc)) 


5ma Siric 


Cac = --{PJof [{naSaUcSc) " {UaSa) iUcSc)) 


Sqab Spc 


D(ab)c = —P'^ K {{naSaUtSbnc) - {UaSaribSb) iUc)) 


Sqab Sruc 


E(ab)c = -P^Jo {{naSanbSbUcSc) - {UaSaUbSb) {UcSc}) 


Spa Srric 


Cac = -P'^KJo {{naUcSc) " ijla) {UcSc)) 



TABLE II: Difi^erent contributions to the fluctuation functional of the Random BEGC model in the replica formalism. All the 
averages have to be evaluated at the saddle point given by Eqs. The subscript {ab) means distinct pairs of indexes. 

I 

A.l The Stability of the RS Solution the stability of the RS solution. The RS Ansatzis imple- 

mented by the substitutions: 

Inserting in the above expressions the Ansatz that all 
replicas are equivalent (same density, same magnetiza- 
tion and same overlap between all of them) we can study qab = (1 ~ Sab) Qo ] Pa = P ] — fn ■ (124) 



+ ^ Spa Da{cd) Sqcd + ^ Spa Bac Spc 



a^{cd) 



Spa Fac Snic + ^ Snia Ea{cd) Sqcd 



a,{cd) 



+ ^ Snia Fac Spc + ^ Sma Cac Snic 



and the eigenvalues equations, then, come out to be 
Sqab + ^ A Sqcd (121) 

(a6),(cd) 

+ X] ^{ab)c Spc + ^ E(ab)c Smc = A Sqab , 



PK Spa + D<cd) Sqcd (122) 
c,d 

+ ^ Bac Spc + ^ Fac Suic = A Spa , 

c c 

PJo Sma + ^ Ea^cd) Sqcd (123) 

c,d 

+ ^ Fac Spc + ^ Cac Snic = A Snia ■ 
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rtcpilCd oynilllCTIlC Jl;XpicoolOIlo lOI XIlc V^OcIIICIcIKd 111 




A{ab)(ab) = 


-I3^[{nant} - (naSaUtSbf] = ~P^[J Vy ~p\y) - (/ Vy rh\y)f] 




A(ab)(ad) ~ 


-/3'*[(n,n656n,S',> - {naSaU^Sb) {uaSaU^Sc)] = Kv) "^'(y) ^ (/ m\v)f] 


= Ai 


A{ab)(cd) = 


-(3''[{naSanbSbncS.niSi) - {uaSaUbSb) {ucSrUiSd)] = m^(y) - (/ Vy m\y)f] 


= Ao 


Baa = 


-P^k^[{na) - {na) {ub)] = - {jVy ~p{y)f] 


= Bi 


Bac ~ 


-P^K^ [{nan,) ~ {na) (n,)] = -/J^if^ [/ Vy p\y) - (/ Vy 


= Bo 


Caa = 


-/3Vo2[K) - (n„S„)'] = ~l3\Jl[j Vy p{y) - (/ Vy m(y))'] 


= Ci 


Cac = 


-/JVo^ [KS„7i,5e) - {naSa) {ucS,)] = p\y) - (/ m(y))'] 


= Co 


D(ab)a = 


-(3^k[{naSanbSb) - (n.SanbSb) (rz,>] = ~P^k[^ Vy n?{y) - J Vy rh''{y) J Vy p{y)] 


= Di 


D(ab)c ~ 


-f3''k[{naSantSbn,) - {uaSan^Sb) (n,)] = ~l3^k[J Vy p{y) m\y) - j Vy rh\y) J Vy p{y)] 


= Do 


E{ab)a ~ 


~P'^Jo[{nanbSb) - {uaSanbSb) {uaSa}] = -/9^Jo[/ Vy p{y) m{y) - J Vy m'iy) J Vy m{y)] 


= Ei 


E(ab)c ~ 


-l3^Jo[{nanbSbncSc) - {uaSanbSb) {n^Sc)] = -/3Vo[/ Vy m^{y) - J Vy m^{y) J Vy m{y)] 


^Eo 


Faa — 


^l3^Jok[{naSa) - {na) {naSa)] = -~P^Jok[J Vy m{y) - j Vy p{y) J Vy m(y)] 


= Fi 


P — 

^ ac — 


-P^Jok[{nan,S.) - {na) {n^S,)] = -P^Jok[J Vy p{y) m(y) - J Vy p{y) J Vy m{y)] 





TABLE III: Here we report the coefficients of the second order term in the expansion of the free energy functional lllUII of the 
Random BEGC model around the RS solution. The left hand side of equalities shows the generic expression of the coefficients, 
whereas the right hand side gives the coefficients evaluated with the Ansatz of replica symmetry. On the far right column the 
abbreviation for the RS coefficients are listed. 



Using the definition 



Go = l3Kp + 



P 90 



(125) 



the replica one site Hamiltonian pi2|) becomes: 

-/?i/'(p,m,go) = eo^na (126) 



+/3 ( Jom + M X! 



With this measure we have to compute the averages 
(oqj . . . Oa^) occurring in the saddle point equations pU|) - 
(|13|l and in the coefficients of the eigenvalues equations 
(see Tab. HD. 

"Opening" , with the Hubbard-Stratonovic method, 
the squared sum in Eq. (|126() into the Gaussian inte- 
gral of an exponential with linear exponent we get 



-I3H' _ 



Vy exp 



with the one-index Hamiltonian 

Ha{y) = QoUa + h{y) UaSa 



(127) 
(128) 



and 



Hy) = y — — ^ Jom + h 



(129) 



We also define the one-index partition sum as 

Zaiy) = J2 e-^"-^^y^ = 2 + 2 coshh{y) , (130) 

and the one-index average 



Zc{y) 



(131) 



that is a function of y. 

With the help of Eqs. (fnHji . ifT^ we are now able to 
compute the nominator of the average {oai ■ ■ ■ Oa^)- 



-f3H' 



{S}.,{n} 



/oo k, n 

n E o^-'"-''' n 



(132) 



/OO 
-Dy [Zi{yrl[(^'^-)a. , (133) 



1=1 Sa, -n, _ 



18 



so that the complete expression of the average is eventu- 
ally given by: 

lim (oai ... Oa J = 1™:^ / 'Dy [Zi{y)]"T\{oa^) 

/oo 
n ' (134) 

-°° i=l 

since Z' = J^^Vy [Zi{y)]^ ^ 1 in the zero replicas 
limit. 

In our case Oa can be any combination of Ua and UaSa 
occurring in Eqs. 1114() - (|116() and in the coefhcients of 
the second order expansion term of the thermodynamic 
potential [Eq. ((TTTjl ] reported in Tab. HH In Eq. ifT^ 
Ua is coupled to Qg, whereas UaSa is coupled to h. Their 
one-index averages can, then, be obtained as: 



dlog Za{y) coshh{y) _ 
(n„)„ = — = — ^— = p(y) 



960 e^^o + cosh/i(2/) 



(135) 



, c \ 51og^a(2/) sinh/i(y) ^ 



dh e~®o + coshh{y) 



(136) 



Using these results we get Eqs. ()19|I - H21|) and the ex- 
pression of the coefficients of S^G for the RS solution. 
We report the complete list in table IIIll Filling in those 
terms, the sum occurring in the eigenvalues Eqs. (|121(l - 
(|123|l can, thus, be written in the RS scheme as: 

J2 Aab){cd) Sqcd = {A2 - 2Ai + Ao) 6qab (137) 

(cd) 

+ {Ai - AQ)''^{bqac + ^<?fcc) + X! ^'i'"'- 
c cd 

X! ^{ab)c ^Pc = {Dl - Da){Spa + 5 Rb) -|- i^O ^P'' 



(138) 



E{ab)c Srric = {El - Eo){5ma + Srub) + Eq^St. 



(139) 



X! ^Q(cd) Sled = {Dl - Do) ^ Sqa 



Dn 



J2Sqc 



(cd) 



cd 



(140) 

^ Bac 5pc - (Bi - Bo)5pa + BoY^ (141) 

c c 

^ Fac Smc = {Fi - Fo)Sma + FoJ^ (142) 



X! ^a(crf) ^Icd = {El - Eo) ^ 6qac + X! ^'^'"^ 



(cd) 



cd 



Fac 5pc = {Fi - Fo)6pa +FoY, 



(143) 
(144) 



Cac Srric = {Ci - Co)5ma + Cpy^Sr 



(145) 



In the case Jo = (no ferromagnetic phase) Cab = 
E[ab)c = Fab = 0, Va, 6, c and we are left with a system 
of equations for n{n— l)/2 variables 5qab and n variables 
5pa {Snia are not involved, since always coupled with 
Jo). In total, the dimension of the space of solutions is 
rftot = n{n + l)/2. 

• Ao 

To obtain the first eigenvalue Ao we analyze the 
equations l|121() - (|122|l in the subspace 



Sqab = , Vfe 



Spa^O ..yb . (146) 



These are 2n equations and the dimension of this 
subspace is = ^^tot — 2n = n(n — 3)/2, corre- 
sponding to the degeneracy of Ao. Eq. Ij 146(1 also 
imphes that X^afc ^l^b = J2a ^Pa = 0- 

Using the conditions (|146|l . the equations (|121|l . 
p22|l are easily reduced, in this case, to the sin- 
gle equation 

{P^ + A2 ~ 2Ai+Ao) Sqab = kii Sqab, (147) 

yielding 

Ao = -f A2 - 2Ai + ^0 . (148) 



• Ai 



We look at the second eigenvalue in the subspace 
Y^Sqab^Q ; ^5pa = 0, (149) 

ab a 

where, as opposed to the previous case, ^lab ^ 
0, V6 and SaPa 7^ 0, Va. The Ao and the Ai sub- 
spaces are orthogonal. 

Subtracting two [the equalities given in Eq. (|149|l ] 
and the degeneracy of the Ao to the total dimension 
we get the degeneracy of Ai : = dtot — dAo — 2 — 
2{n~ 1). 

Eqs. (|121|I - H122() reduce in this subspace to 

[(3^ + A2 + {n- A)Ai -{n- 3)Ao] J2 ^lab (150) 

6 

+ {n - 2){Di - Do) Spa = Ai ^ Sqab 

b 

{Di~Do)Y^1ab (151) 

6 

+ (/3X + Bi- Bo) Spa = Ai Spa 

for which the eigenvalue comes out to be 
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Ai = - 
^ 2 



P'^+A2 + {n- 4)Ai -{n~ 3)Ao + fiK + Bi - 
1 



(152) 



±- 



/32 + ^2 + - 4)Ai - (n - 3)ylo - {I3K + Bi - Bo) + 4(n - 2){Di - D^f 

I 



• A2 

To find A2 we are left witli tlie subspace of solutions 
such that 

^ 5qab ^0 ; ^ ^ . (153) 



and whose dimension is 2. 

The eigenvalues equations become now 



0'+A2 + [n - 2) ( 2Ai + ^^^0 



ah 



[n - 1) [2D^ + {n- 2)Do] ^ Spa - A2 ^lab 

a ah 

(154) 

Di + '^^D^Y.^qab (155) 

^ ab 

f3K + Bi + {n- l)Bo) ^ 5pa = A2 ^ Spa 



so that it, finally, holds 



13^ + A2 + 2{n - 2) Ax + 



(n- 2)(n- 3) 



n-2 



-Do 



(156) 



±- 



/32 + A2 + {n - 2) ( 2Ai + 



n-2 



-Do 



+ 2{n-l)[2Di + in-2)Do] 



In the limit for n ^ the eigenvalues Ai (|152|l and A2 p56|l are degenerate. They reduce both to 



Ai A2 



A2 - AAi + 3Ao + f3K + Bi-Bo 



2 



f3^+A2- AAi + 3^0 - PK-Bi 
I 



-8(i?i-i?o)2 



(157) 



If we also put /i = we simplify things much, since in 
this case the stable solution is qo = 0. This brings to 



m{y) = , 



Po , 



(158) 
(159) 



and, therefore, no integral in the coefficients of table Hill 
has to be carried out anymore. This leads to Ai = Ao — 
Di ^ Do = Bo = and 



A2 
Bi 



-/3'i^Vo(l - Po) ■ 



(160) 
(161) 



obtain 



ko = P\l-l3'pi) , 



Ai 



I3K 



l-fikpo{l~Po) 



(162) 
(163) 



The above analysis is valid for K > Q. This is always 
the case ik the biquadratic coupling \s K > Q. When a 
negative K occurs, however, at high enough temperature 
(T > l/(2|_fi'|)) K becomes negative. Taking this into 
account and repeating the whole scheme of computation 
the final result for if < is, in the case of out interest, 



Ai 



(3k l-l3kpo{l-po) 



(164) 



Substituting these coefficients into Eqs. p48|) - (|157|l we 
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APPENDIX B: Small q Expansion 

The expansion to the fourth order in r — p — po and 
Qab is: 



(165) 



Pfi{qab},r)=Pfo-PK[po-j)r 



— po(l-po)(l-2po)?- 



A 
24 



Po(l-7po + 12p'-6pV 



Ao 1 , (djy ~ ^, 
+ -^-^^Kplil- p^) 



1, 



-Tr q2 



48 



\l + 3PK{l~ po)r\ -Tt q3 
L J n 



X! 9afc + 6P0 X! ^^fc 

(a6) (a6c) 



lac 



2 X ^ 2 2 2 \ ^ 



b Qbc (led Qda 



{abed) 



{abed) 



where ^ is a sum over distinct indexes. The para- 
magnetic contribution is 



Nb + I 

i=0 
Nb 

] 1^ [ ^ab - ^ab 



b 



(169) 



(170) 



with qNB+i = 9-1 = 0, rriNB+i = 1, mo = n, e^^ = 5ab, 
ef^ a, n X n matrix with all elements equal to 1. 

The matrix e'"^^ has mi/rrii^i blocks along the diagonal 
and its summation rules are: 



/_^'^abHe — "'-niax(ij)tac 
6 

V 5 J^'> - e(^') 



E4' 



b - 



b^be ^cd 



min(i,j ,k) 
mmax(i,j,k) ''TT'2°'^nistx{i,i,k)^ad 



E4' 



b ^be ^cd ^de 



(171) 
(172) 
(173) 
(174) 

(175) 



bed 



™ ^min(i,j,k,l) 

— "lmax(i,j,k,l) "l2"<ij„ax(i,j ,k,l) "^3'-'imax(i,j,k, 1)^06 

where ^ is a sum over all indexes. 



Term Trq^ = J^ab 1'^'> ^^<^ 



/3/o = ^P^ + log(l-po) 



and 



1 + cxp 



-Pp - ISKT/J 



(166) 



(167) 



is the value of the paramagnetic density evaluated along 
the second order phase transition line. 

The above expansion is valid in the neighborhood of 
the second order transition line, Aq = 0, down to, and 
including the tricritical point (Aq = U Ai = 0). 

The various terms appearing in the expansion of the 
free energy functional, Eq. (|165() are expressed in the 
following. In order to compute it in a generic A'^g-RSB 
scheme we just need the following expression for the form 
of the overlap matrix qab 



i=NB 



Qab 



-qWB^ab + E ^ '^'^-l) ^ib + 50 (1^^) 



Nb 

E 9af> Ibe = 
b 1=0 

Nb 

E 9afe Iba = E ~ ^i+l) 

b j=0 



Y.q!{m,el^~m,+,el:+'^) ,(176) 

(177) 



In the full replica symmetry breaking limit and in the 
zero replicas limit the breaking parameters become con- 
tinuous between and 1: 



rrii i ruk — rrii^i — > —dx 



(178) 



Moreover the same structure of Eq. p77|) holds for any 

Eah 9a6' thus 



(179) 
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type of sum 


contribution 


multiplicity 


A: i < j < k 





6 


B: i = j < k 


qk ql{mk - mk+i)imi - rut+i) 


3 


C: i < j = k 





3 


D: i = j — k 


ql \{mk - TUk+i)^ - mk+i{mk - mk+i) 


1 



TABLE IV: Contributions to the trace of q for a generic number of breakings of the replica symmetry: only two kind of terms 
are different from zero. 



Term Trq^ = J^abc 1<^i> 9=' 



Using p7()|l this trace can be written as 

E 



be 



'^abHc'^ca '^ab He '^ca '^ab'^be '^ca ^'^ab He '^ca 



(180) 



"'ob He "^ca 



ab be ca 



^abHc ^ca 



^ab He ^ca 



In table llVl we report the four possible non degenerate contributions coming out from this sum and their multiplicity. 
In the full replica symmetry breaking limit (and for number of replicas n Q): 

rui ^ y mi - m^+i -dy . 
Thus the above expressions for the terms contributing to the sum reduce to 

1 px 

dx q{x) I dy q{yf , 



B 
D 



dx X q{x) 



where we have neglected terms of order {dxY . Summing up, the trace comes out to be: 



(181) 
(182) 



lim -Tr = 2,B + D 

n^O n 



3 dx q[x) I dy q{y)'^ + dx x q{x)^ 
Jo Jo Jo 



(183) 



Term Trq'' = J2abcd H": <l<^<i 1<i<i. 
Exploiting the formulation of Eq. 1)190(1 once again we can obtain 



0, Nb 

^ Qab qbe qed Qda ^ ^ li 1] Ik ^ {^ab He "ed "da 

1. j,k bed 



bed 



(i) _ (i+i)(j)(k) (i) _ A^)(j+i)(k) (i) ^(i),(j),(fe+i),(0 

fc„i,t;,„ t„j fc^„ fc^h '^bc'^ed'^da '^ab'^bc '^ed'^da 



He '^cd 



"da 



Ji) Aj) Ak) (1+1) , Ai+1) Aj+1) Ak) M) , Ai'+i) Aj) Ak+1) M) , Ai+i) Aj) M M+i) 

"ab "be "ed "da '^ab "be "cd "da ~^ "ab "be '^ed '^da ^ "ab He '^ed "da 



+H,b<^ 



{i)Aj+i)Ak+i)Ai) : Ai)Aj+i)Ak)Ai+i) . A^) Aj) (k+i) AW) (i+i) Aj+i) (k+i) M) 



be "ed "da ~^ ^ab He 



^cd ^da 



^abHc ^ed 



"da 



"be "cd 



"da 



Ai+1) Aj+1) Ak) (i+i) _ Ai+1) Aj) Ak+1) (i+i) _ A':)Aj+i)Ak+i)Ai+i) , A''+i) Ai+i) (k+i) AW) 

"ab "be "ed ''da ''ab ''be ''ed ''da ''ab''bc ''ed ''da ^"-ab He ''ed ''da 



Contribution and multiplicity are reported in tabled For Nb — * oo and n ^ the terms are 

B^ - I dx q{x) I dy q{y) r dzq^{z) (184) 
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type of sum 


contribution 


multiplicity 


A: i < j < k < I 





24 


B: i = j < k < I 


qi qk qf{mi — mi+i){mk — mk+i){mi ~ rm+i) 


12 


C: i < j = k < I 





12 


T): i < j < k = I 





12 


E: i= j <k = l 


q^qlinik - mfe+i)^(mi - mi+i) 


6 


F: i = j = k < I 


qi qi{mi ~ rrii+i) [(mfc - m^+i)^ - mk+i{mk - ruk+i) 


4 


G: i < j ^ k = I 





4 


R:i = j = k^l 


qf{mi - nii+i) \m'^^-^ - mi+i{mi - m^+i) + {rrii - rrii+i)^ 


1 



TABLE V: Contributions and multiplicities of the terms in the sum of Eq. H184^ . Only for kinds of terms are non-zero. Of 
them, term E turns to be of order (dx)^ once the zero replica limit has been performed. 



F - 
H 



dx q{x) / dy y q (y) 



dx X q (x) 



(185) 
(186) 

(187) 



Term Y.abc'&bllc 
Eventually this last object reads: 

51 '^^f-'^^c ^^llbY^ lie 

be be 
i 3 

the FRSB hmit (for n ^ 0) being 



(188) 



abe 
1 



dx q^{x) I dy q^{y) 



dx q (x) 



1. Full RSB free energy for disordered BEG around Tc 

Computing the expansion for small q and p ~ poi Eq, 
p65|l . in the FRSB scheme, yields 



+ / dx X q{xY 



{f3Jf 



-Po 12 / dx q{x) / dy q{y) / dz q{z) 



+4 / dx q{x) I dy y q{y)-^ + / dx x^q^{x) 



+ ^^fF~/'o (1 - 3po)^ / dx q{xf 



{(iJf 



■Po (1 - 3po) / dx q\x) 



1440 



-pI{1- 30po + 285pg - 900pg 

+900p^) / dx q^{x) 



Ao 



dx q{xY 







The expansion is valid along the critical line p = T 
down to the tricritical point {pc,Tc) given by Eqs. I|29t - 
(189) EOJ. 

The last term is introduced for the case K = 0, i.e. 
for the Ghatak-Sherrington model, in which 1 — 3po goes 
to zero at the tricritical point p^ = 1/3 (e.g., decreasing 
temperature or density along the second order line). In 

this case the relevant quartic term Jq dx q{x)^, respon- 
sible for the replica symmetry breaking, vanishes. This 
means that, in this case, the symmetry breaking is weaker 
than, for instance, in the SK case. 

Putting p ^ 1 and r = in the above expression we 
obtain the same result of Ref. 'Hl'l , where the same kind 
of expansion for SK model is performed, up to the forth 
(190)order. To allow a straightforward check we rewrite the 
above formula in the way of Ref. 



my 

24 



{(3Jf 



Po(l - Po)(l - 6po + 6po)r"* 



/?/ = /9/o 

-M!(l-p^(/3J)^) fdxq{xf 



(191) 



/3Kp2(l-po) 



— (2- 3po)r 



dx q{xy 



Po / da; a; q{x) + 3q{x) / dy q{y) 



1 + 3f3k{l - po)r 



dx q{x) / dy q{y) + 



24 



8 rl 



dx 



^pU^ + 2-00(1 - 3po)' 



q{xY 
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-2p^(l-3po) q{xY / dyq{yY 
Jo 

+ 12pt q{x) I dyy q{yf 
Jo 

rx py 
+36po q{x) / dy q{y) / dz q{zf 
Jo Jo 

+0{q'^) + 0{r) 

Variation with respect to q{x) and further differentia- 
tion with respect to x lead to an integral equation that 
can be reduced to 



Qxp 



q{x,p) X {\ — 2ipY -\- Qxp"^ ^ ^ 



d^q{x,p) = (193) 
Below the critical temperature the solution is: 

Ci X 



q{x,p) = 



(1 - 3/9)2 ^(1 _ 3p)2 + Qp2^2 



(194) 



as far as a; < Xm{p), otherwise q(x,p) = qi for x > 
xm[p)- The value xm is given by equating qi = q{xM , p)- 



XM^p) 



gi(l-3p)3 



gi(l-3p)3 



^Cf-6p^l-pYql 



Oiq') 
(195) 

Notice that in Eqs. ifTO^ . ifT^ p po + and the 
formulas have yet to be expanded in r. The above com- 
putation can be simplified neglecting the quartic terms in 
Eq. (|165() which are irrelevant with respect to the RSB, 
i.e. all but the one involving g^^. In this case q{x) is 
simply 



q{x) — A{r) X 



(196) 



with 



for X < Xm ^ qi/A{r) and q{x) — qi for x > Xm- If 
r and po ^ 1 ^(0) reduces to 1/2 (SK model),^^ 
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